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ABSTRACT 

We present the first coordinated VLBA/VLTI measurements of the stellar diameter 
and circumstellar atmosphere of a Mira variable star. Observations of the v = 1, J = 
1-0 (43.1 GHz) and v = 2, J = 1 - (42.8 GHz) SiO maser emission toward the 
Mira variable S Ori were conducted using the VLBA. Coordinated near-infrared K- 
band measurements of the stellar diameter were performed using VLTI/VINCI closely 
spaced in time to the VLBA observations. Analysis of the SiO maser data recorded 
at a visual variability phase 0.73 show the average distance of the masers from the 
center of the distribution to be 9.4 mas for the v = 1 masers and 8.8 mas for the 
v = 2 masers. The velocity structure of the SiO masers appears to be random with 
no significant indication of global expansion/infall or rotation. The determined near- 
infrared, iT-band, uniform disk (UD) diameters decreased from ~ 10.5 mas at phase 0.80 
to ~10.2mas at phase 0.95. For the epoch of our VLBA measurements, an extrapolated 
UD diameter of 9{f D = 10.8 ± 0.3 mas was obtained, corresponding to a linear radius 
of R* D = 2.3 ± 0.5 AU or i$ D = 490 ± 115 R Q . Our coordinated VLBA/VLTI 
measurements show that the masers lie relatively close to the stellar photosphere at 
a distance of ~ 2 photospheric radii, consistent with model estimates. This result is 
virtually free of the usual uncertainty inherent in comparing observations of variable 
stars widely separated in time and stellar phase. 

Subject headings: masers — stars: AGB and post AGB — stars: atmospheres — stars: 
late type — stars: mass-loss — techniques: interferometric 
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1. INTRODUCTION 

The evolution of cool luminous stars, including Mira variables, is accompanied by significant 
mass loss to the circumstellar envelope (CSE) with mass- loss rates of up to 10~ 4 M /year (e.g. Jura 
& Kleinmann 1990). There are currently a number of tools that can be used to study the stellar 
surface and the CSE at various wavelengths. Optical/near-infrared long-baseline interferometry has 
provided information regarding the stellar diameter, effective temperature, center-to-limb intensity 
variation, and the dependence of these parameters on wavelength and variability phase for a number 
of Mira variables (e.g., Haniff, Scholz Sz Tuthill 1995; van Belle et al. 1996; Perrin et al. 1999; Young 
et al. 2000; Thompson, Creech-Eakman & Akeson 2002a; Thompson, Creech-Eakman & van Belle 
2002b; Woodruff et al. 2004). The structure and dynamics of the CSE of Mira variables and 
supergiants have been investigated by mapping the SiO maser emission at typical distances of 2-4 
stellar radii toward these stars using very long baseline interferometry (VLBI) at radio wavelengths 
(e.g., Boboltz, Diamond & Kemball 1997; Kemball & Diamond 1997; Boboltz & Marvel 2000; 
Hollis et al. 2001; Diamond & Kemball 2003). Dust shells outward of the SiO maser emission 
zone have been studied using mid-infrared interferometry (e.g., Danchi et al. 1994; Greenhill et 
al. 1995; Townes 2003). Near- and mid-infrared interferometry have also been used to study the 
warm extended atmosphere of Mira variables containing gaseous H2O, SiO and CO molecules (e.g., 
Perrin et al. 1999; Mennesson et al. 2002; Ohnaka et al. 2004). 

Theoretical models describe and predict photospheric center-to-limb intensity variations in- 
cluding the effect of close molecular layers (e.g., Bessel, Scholz &: Wood 1996; Hofmann, Scholz & 
Wood 1998; Tej et al. 2003; Ireland, Scholz & Wood 2004), the location and dynamics of SiO maser 
emission zones (Humphreys et al. 1996, 2002), as well as gas and dust components of the stellar 
winds (e.g., Winters et al. 2000; Hoefner et al. 2003). 

Multi-wavelength studies of the stellar surface and the CSE, utilizing a combination of the 
above techniques, are well suited to study the structure of the CSE, the mass-loss process, and 
the evolution of these stars. It is common for infrared observations of the stellar photosphere 
or the circumstellar dust to be compared to radio observations of circumstellar masers from the 
literature and vice versa. For example, Danchi et al. (1994) compared the inner radii and extents 
of dust shells as determined from Infrared Spatial Interferometer (ISI) observations with estimates 
of the photospheric radii and measurements of the extents of the SiO, H2O, and OH maser shells. 
Similarly, Monnier et al. (2004) compared high-resolution data on the stellar diameters and dust 
shells of two stars (VX Sgr and NML Cyg) from Keck aperture masking and Infrared Optical 
Telescope Array (IOTA) interferometry with previously published maps of the SiO maser emission 
toward these stars. Such comparisons, however, are somewhat limited due to the inherent variability 
in the CSE and the star itself. The SiO masers for example have been shown to have proper 
motions and intensity variations on time scales of a few weeks (Boboltz, Diamond & Kemball 1997; 
Diamond & Kemball 2003) . Similarly, interferometric measurements of the stellar parameters vary 
with stellar variability phase as discussed above. This inherent variability of Mira stars and their 
surroundings necessitates a more coordinated approach to multi- wavelength studies of these objects. 
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The first such coordinated multi-wavelength study of a late-type variable was performed by 
Greenhill et al. (1995) for the star VX Sgr. Coordinated measurements included both mid-infrared 
observations to determine the extent of the circumstellar dust shell and VLBI mapping of the SiO 
maser zone. More recently, Cotton et al. (2004) compared Very Long Baseline Array (VLBA) 
images of the SiO maser emission of several Mira stars with near-infrared diameters obtained at 
the Infrared Optical Telescope Array (IOTA) interferometer (Mennesson et al. 2002), and with 
literature values of the mid-infrared dust shell extensions. 

Here, we present the first results from a series of coordinated VLTI/VLBA observations of late- 
type variable stars, in particular, for the Mira variable S Ori. Our results include the first radio 
interferometric images of the SiO masers (VLBA) and the most accurate photospheric diameter 
measurements using near-infrared i^T-band interferometry (VLTI/VIN CI) . 

2. CHARACTERISTICS OF S ORI 

The target of our observations, S Ori, is a Mira variable star with spectral type M6.5-M9.5e 
(Kholopov et al. 1985) and V magnitude 8.4-13.3 (Pojmanski 2002). Whitelock, Marang & Feast 
(2000) report a quite low near-infrared J-band pulsation amplitude of 0.55 mag. They also derive 
a mean bolometric magnitude m\> \ = 3.08, with an amplitude Am-boi = 0.49. Over the last 100 
years, the variability period of S Ori has been seen to randomly vary between about 400 and 450 
days (Bedding et al. 1999; Whitelock, Marang & Feast 2000; Merchan Benitez & Jurado Vargas 
2002). For our study, we adopt a stellar period of 420 days and a time of last maximum visual 
brightness of JD 2452331, as given by Pojmanski (2002). This is consistent within ± ~ 5 days with 
an inspection of the most recent AAVSO and AFOEV data, van Belle et al. (1996) measured the 
near-infrared .KT-band uniform disk (UD) diameter of S Ori at variability phase 0.56 to be 10.54 ± 
0.68 mas, and derived a Rosseland mean diameter of 11.70 =L 0.75 mas. We adopt the distance to 

5 Ori to be 422 ± 37 pc, as given by van Belle et al. (1996), based on measurements by Wyatt 

6 Cahn (1983) and Young (1995). This is consistent with the value of 430 pc from Whitelock, 
Marang & Feast (2000) based on their K period-luminosity relation. S Ori exhibits SiO and OH 
masers (Benson et al. 1990), but no detection of any H2O maser emission has been reported. Sloan 
& Price (1998) report a relatively low dust emission coefficient {DEC), i.e. the total emission of 
the dust to the total emission of the star in the wavelength range from 7.7-14.0 /im, of 0.24 for S 
Ori (for comparison, DEC = 0.46 for oCet). 
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3. OBSERVATIONS AND REDUCTION 

3.1. Near-infrared K-band interferometry 

We obtained near-infrared .fT-band interferometric measurements of S Ori using the ESO Very 
Large Telescope Interferometer (VLTI) equipped with the commissioning instrument VINCI and 
the two VLTI test siderostats (effective apertures ~40cm). VLTI stations D 1 and B3, forming an 
unprojected baseline length of 24 m, were used for 13 nights between 2003 January 25 (JD 2452665) 
and 2003 March 31 (JD 2452730). The S Ori visual variability phase ranged from 0.80 to 0.95 for 
these dates. Available data taken in 2002 December using an 8 m ground baseline and in 2003 
September using a 16 m ground baseline were not used because S Ori was only marginally resolved 
with these shorter baselines (|V| > 0.9). A recent general description of the VLTI can be found 
in Glindemann et al. (2003) and references therein. During each observation night, several series 
of typically 500 interferogram scans were recorded on S Ori as well as on several calibration stars. 
The calibration stars used and their adopted properties are listed in Table 1. For S Ori, an effective 
temperature of 2500 K was used, which is consistent with the value of 2313 ± 110 K obtained by 
van Belle et al. (1996) for phase 0.56. We note that variations of the effective temperatures up to ± 
500 K result in variations of the squared visibility amplitude of less than 0.3%. The scan length for 
the S Ori and all calibration stars' observations was 230 //m, the scan speed was 616//m/sec, and 
the fringe frequency, corresponding to the time to scan one interferometric fringe, was 289.5 Hz. 

Mean coherence factors were obtained for each series of interferograms using the VINCI data 
reduction software, version 3.0, as described by Kervella et al. (2004), employing the results based 
on wavelets transforms. Calibrated squared visibility values for S Ori were obtained by calibration 
of the mean coherence factors as described in Wittkowski et al. (2004), with a time kernel of 5 hours 
to convolve the measured transfer function. The errors of the calibrated S Ori squared visibility 
values include the scatter of the single-scan coherence factors, the adopted errors of the diameters 
of the calibration stars, and the observed variation of the transfer function during each night. The 
calibrated visibility values are listed in Table 2 and are also available in electronic form from the 
authors upon request. Table 2 also lists the date and time of observation, the spatial frequency, the 
position angle of the projected baseline (East of North), the calibrated squared visibility value V 2 
and its error ay2, as well as the number of processed interferograms. The effective wavelength is 
2.19/im. All observations were performed at similar angles of the projected baseline ranging from 
72° — 74° E of N. The squared visibility values for S Ori were grouped into four bins of five nights 
each (bin width ~ 1% of the visual variability period of 420 days). 

We characterize our measurements with a best-fitting uniform disk (UD) diameter for each of 
our four epochs. Our uniform disk model takes into account the broad passband (1.9-2.5 ^m) of the 
VINCI instrument, as described in Wittkowski et al. (2004). Figure 1 shows the obtained squared 
visibility amplitudes of S Ori together with the best fitting uniform disk models for each of our 
four VLTI/ VINCI epochs. The larger scatter of the visibility data for the March 16-21 epoch, as 
compared to the other epochs, can be explained by a larger time difference of the calibration stars. 



- 5 - 



Table 3 lists the observation dates and the determined K-band UD diameters and their associated 
formal errors. Additional calibration uncertainties are estimated to be ~ 0.1 mas, based on an 
analysis of sub-volumes of the data and differing methods of interpolating the transfer function. 
The addition of a possible circumstellar dust shell is not expected to affect the K-band UD diameter 
significantly (see for instance discussions in Woodruff et al. 2004; Ohnaka et al. 2004). In particular, 
Ohnaka et al. (2004) modeled an optically thin dust disk for a similar Mira star, RR Sco (dust 
emission coefficient 0.21, compared to the value of 0.24 for S Ori), for the purpose of comparison to 
mid-infrared interferometric data. At a wavelength of 2.2 /xm the scattered and thermal emission 
from this dust disk is lower by a factor of ~ 1000 compared to the attenuated star light and the 
effect on the squared visibility values at V 2 ~0.5 is less than 0.2% (Ohnaka, K. 2004, private 
communication). This would lead to an effect on the uniform disk diameters in Tab. 3 of less than 
0.02 mas. Discussed in Sect. 4.1 below is the relationship between a K-band UD diameter and a 
physically more meaningful angular size of the stellar photosphere. 



3.2. SiO Maser Observations 

We simultaneously observed the v = 1, J = 1 - (43.1 GHz) and v = 2,J = l — (42.8 
GHz) SiO maser transitions toward S Ori (a = 05 A 29 m 00 s .9, 5 = -04°41'32".7, J2000) on 2002 
December 29for 6 hrs starting at 02:43 UT (JD 2452636.6, visual variability phase 0.73). S Ori 
and a continuum calibrator (0359+509) were observed using the 10 stations of the VLBA. The 
VLB A is operated by the National Radio Astronomy Observatory (NRAO). 1 Reference frequencies 
of 43.122080 and 42.820587 GHz were used for the v = 1 and v = 2 SiO transitions respectively. 
Data were recorded in dual circular polarization using two 8-MHz (56.1 km s _1 ) bands centered on 
the local standard of rest (LSR) velocity of 18.0 km s -1 . System temperatures and point source 
sensitivities were on the order of ~150 K and ~11 Jy K _1 respectively. 

The data were correlated at the VLBA correlator operated by NRAO in Socorro, New Mexico. 
Auto and cross-correlation spectra consisting of 256 channels with channel spacings of 31.25 kHz 
(-0.2 km s" 1 ) were produced by the correlator. Calibration was performed using the Astronomical 
Image Processing System (AIPS) maintained by NRAO. The total intensity data were calibrated 
in accordance with the procedures outlined in Diamond (1989). The bandpass response was de- 
termined from scans on the continuum calibrator and used to correct the target source data. The 
time-dependent gains of all antennas relative to a reference antenna were determined by fitting 
a total-power spectrum (from the reference antenna with the target source at a high elevation) 
to the total power spectrum of each antenna. The absolute flux density scale was established by 
scaling these gains by the system temperature and gain of the reference antenna. Errors in the gain 
and pointing of the reference antenna and the atmospheric opacity contribute to the error in the 



x The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 
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absolute amplitude calibration, which is accurate to about 15-20%. 

To correct any instrumental delay, a fringe fit was performed on the continuum calibrator 
scans, and residual group delays for each antenna were determined. Variations in the residual delays 
ranged from 2-4 ns resulting in phase errors of no more than 1.5-3° across the 8-MHz band. Residual 
fringe-rates were obtained by fringe-fitting a strong reference feature in the spectrum of each maser 
transition. For both transitions we used the same channel at a velocity Vlsr = 16.2 km s _1 . 
The resulting fringe-rate solutions were applied to all channels in each spectrum respectively. An 
iterative self-calibration and imaging procedure was then performed to map this reference channel 
for each transition. The resulting residual phase and amplitude corrections from the reference 
channels at 42.8 and 43.1 GHz were applied to all channels in the respective bands. 

In order to accurately compare the distributions of the two maser transitions, it is desirable 
to determine a common spatial reference point. However, after the fringe-fitting step to determine 
residual fringe-rates, all absolute position information is lost for the VLBA data. We accomplished 
the registration of the two transitions by applying the fringe-fit solutions from the v = 2, 42.8- 
GHz transition to the v = 1, 43.1-GHz data, re- mapping the v = 1 transition, and comparing the 
resulting images with those produced from the application of the v = 1 calibration itself. Although 
the images resulting from the application of the 42.8-GHz fringe-fit to the 43.1-GHz data were of 
poorer quality, we were still able to determine an offset between the two sets of images. The offsets 
in right ascension and declination computed from images of two different spectral channels, were 
the same to within 0.03 mas. Subsequent phase self-calibration and imaging were performed for 
the v = 1 SiO data using these computed offsets, resulting in spatially aligned v = 1 and v = 2 
image cubes. This same procedure was applied to the continuum calibrator source, 0359+509, 
and the positions derived agreed to within 0.4 mas, thus providing an estimate of the error in the 
registration of the two maser distributions. 

Final images of the SiO maser emission consisting of 1024 x 1024 pixels (~51 x 51 mas) were 
generated using synthesized beams of 0.58 x 0.17 mas and 0.52 x 0.17 mas for the v = 1 and 
v = 2 transitions respectively. Images were produced for spectral channels from 4.7 km s _1 to 
30.7 km s" 1 forming image cubes of 120 planes. Off-source RMS noise in the images ranged from 
7 mJy to 18 mJy. Figure 2 shows the total intensity contour maps of the v = 1, 43.1-GHz (red) and 
v = 2, 42.8-GHz (blue) SiO maser emission toward S Ori. The contours represent the maximum 
pixel value in each image cube over the LSR velocity range from +10.1 km s _1 to +25.4 km s _1 . 
Subsequent analysis of the image data is described in detail in Sect. 4.2. 
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4. RESULTS AND DISCUSSION 

4.1. The photospheric diameter 

Figure 3 and Tab. 3 show our near-infrared K-band uniform disk (UD) diameters of S Ori 
determined from our four VLTI/ VINCI epochs as a function of observation date and visual vari- 
ability phase, as discussed in Sect. 3.1. The K-band UD diameter decreases almost linearly from 
~ 10.5 mas at visual variability phase 0.80 to ~ 10.2 mas at phase 0.95, i.e. by ~3%. It is not clear 
whether the local diameter minimum of ~ 10.0 mas at phase 0.92 in Tab. 3 and Fig. 3 is caused by 
an additional systematic calibration error as discussed in Sect. 3.1, or by a real minimum of the 
stellar size. The minimum of a Mira K-band UD diameter is predicted to lie at visual variability 
phase ~0.9, according to model predictions by Ireland, Scholz & Wood (2004). Our coordinated 
VLBA measurement at visual variability phase 0.73 occurred about one month before the first 
VLTI epoch at phase 0.8, and is hence close to, but not exactly contemporaneous to the VLTI 
epochs. Both, observations of the diameter change of the Mira star SLac by Thompson, Creech- 
Eakman & van Belle (2002b) as well as Mira star model predictions by Ireland, Scholz & Wood 
(2004) show an almost linear change of the K-band UD diameter between variability phases 0.7 
and 0.9, i.e. in the pre-minimum part of the diameter curve and pre-maximum part of the visual 
light curve. Thus, we can linearly extrapolate our K-band UD diameter data at phases 0.8-0.95 
to the VLBA phase 0.73. For phase 0.73, we determine an extrapolated K-band UD diameter of 
9^ D (VLBA epoch, phase = 0.73) = 10.8±0.3 mas. This value and its error correspond to the mean 
and difference of the two linear extrapolations shown in Fig. 3, which used all four points and only 
the two closest points to the VLBA epoch respectively. This UD diameter corresponds to a linear 
radius i$ D = 2.3 ± 0.5 AU or R^ D = 490 ± 115 R Q with our assumed distance to S Ori. 

A UD intensity profile is often not an ideal representation of the true near-infrared center- 
to-limb intensity variation (CLV) of cool giants in general and of Mira stars in particular. The 
CLV of cool giants has been studied by, e.g. Quirrenbach et al. (1996); Burns et al. (1997); 
Scholz (1998); Hofmann, Scholz & Wood (1998); Wittkowski et al. (2001, 2004); Ireland, Scholz & 
Wood (2004); Woodruff et al. (2004). UD diameters of non- variable giants obtained from visibility 
measurements in the first lobe are usually transformed into physically more meaningful Rosseland 
angular diameters using correction factors determined from atmosphere models in the literature 
(e.g. Claret 2003). For the case of Mira variable stars, few atmosphere models are available. 
We used the hydrodynamic atmosphere models from Hofmann, Scholz & Wood (1998); Tej et 
al. (2003); Ireland, Scholz & Wood (2004), to estimate the relationship between the K-band UD 
diameter and the continuum diameter. While the broad-band UD diameter is affected by various 
molecular bands, the continuum diameter is a better estimate of the real photospheric stellar size. 
These models have been constructed for the prototype Mira stars oCet and RLeo and have been 
compared to observations of several Mira stars by, for instance, Hofmann et al. (2001, 2002); 
Woodruff et al. (2004). The general model results are not expected to be dramatically different 
for other Mira stars such as S Ori (Scholz, M. 2004, private communication). The model predicted 
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difference between the continuum and if-band UD diameters is relatively low in the pre-maximum 
region of the visual variability curve as in the case of our observations. At this phase of 0.73, the 
continuum diameter may be smaller than the K-band UD diameter by about 15% (Ireland, Scholz 
&; Wood 2004) . With this assumption, the continuum photospheric diameter for the epoch of our 
VLBA observation would be 6ph t(VLBA epoch, phase = 0.73) ps 9.2 mas. This angular radius 
corresponds to a photospheric radius of i?phot ~ 420 R or i?phot ~ 1.9 AU with our assumed 
distance to S Ori. With a bolometric flux of mboi = 3.1 ± 0.3 (Whitelock, Marang & Feast 2000), 
our value for 0phot corresponds to an effective temperature T c g ss 2670 K. 

4.2. The v = 1 and v = 2 SiO Maser Emission 

Our total intensity images of the 43.1 GHz and 42.8 GHz SiO maser emission (Figure 2) show a 
typical clumpy distribution of the maser spots within a ring-like structure. For the v = 2, 42.8-GHz 
transition the ring is relatively sparse with nearly all of the masers concentrated to the NW side 
of the shell, and a few features to the SE. The v = 1, 43.1-GHz masers, however, form a more 
typical ring-like structure often seen for other stars with SiO masers. Like the v = 2 masers, the 
v = 1 SiO also has a higher concentration of features on the NW side of the shell. The v = 1 maser 
ring appears to be symmetrical, forming a nearly circular shell. The distribution of the v = 2 SiO 
masers is too sparse to allow us to comment on any symmetry for this transition. 

In order to identify and parameterize maser components, two-dimensional Gaussian functions 
were fit to the emission in each spectral (velocity) channel using the AIPS task SAD. Image quality 
was assessed using the off-source RMS noise and the deepest negative pixel in the image. A cutoff 
flux density was conservatively set to the greater of 8cjrms or the absolute value of the deepest 
negative pixel in the plane. Features with flux densities greater than this cutoff were fit with 
Gaussians to determine component parameters. Errors in the positions in right ascension and 
declination of identified features were computed using the fitted source size divided by twice the 
signal-to-noise ratio (SNR) in the image and ranged from 1 //as for features with high SNR, to 
50 /uas for features with lower SNR. 

Since the ~0.2 km s^ 1 channel spacing is sufficient to resolve the masers spectrally, features 
typically appear in multiple adjacent spectral channels. Positions in right ascension and declination 
and center velocities for the masers were determined using a flux-density-squared weighted average 
for features identified in two or more adjacent channels with a spatial coincidence of 0.2 mas (~2/3 
of the geometric mean of the synthesized beam). The flux assigned to the maser averages was the 
maximum single-channel flux density. The maser components identified using this procedure are 
represented by the circles in Figures 4 and 5. In the figures, point sizes are proportional to the 
logarithm of the fitted flux density which ranged from 0.2 to 4.3 Jy for Figure 4 and 0.1 to 5.3 Jy for 
Figure 5. The line-of-sight (LOS) velocity information for the masers is also represented in Figures 
4 and 5. The top panel of each figure shows the spectrum of SiO maser emission ranging from 10 to 
25 km s _1 color-coded by LOS velocity in increments of 2 km s _1 . The bottom panels of Figures 4 
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and 5 show the spatial distribution of the SiO masers plotted with the same velocity color-coding 
as in the top panel, and with the color of the maser representing its corresponding velocity range 
in the spectrum. Comparing the maser component maps with the total intensity contour maps, we 
see that the identified features accurately represent the emission summed over all velocity channels 
in the image cube. From the component maps there does not appear to be any coherent velocity 
structure indicative of global expansion/infall or rotation. We do note that a group of components 
on the western side of the shell shows a velocity gradient with velocity decreasing with increasing 
distance from the star. Such velocity gradients in the SiO maser emission have been observed 
previously (e.g. Boboltz & Marvel 2000; Hollis et al. 2001). 

To compare the size of the SiO maser distribution to the photospheric diameter of the star as 
measured by the near-infrared i^T-band interferometry we determined the average distance of the 
masers from the center of the distribution. To accomplish this, we first determined the center of 
the distribution by performing a least-squares fit of a circle to the combined v = 1 and v = 2 maser 
component data. This fit produced a common center from which we computed the mean maser 
angular distance rsio an d the standard deviation for each transition independently. The mean 
angular distances from center for the observed SiO masers at v = 1, 43.1-GHz and v = 2, 42.8-GHz 
are rsio = 9.4 and rgio = 8.8 mas respectively. These distances are indicated by a dashed circle and 
are listed in the lower panel of Figures 4 and 5. Standard deviations of the distances are cisio = 1-4 
and cisio = 1-7 mas for the v = 1 and v = 2 transitions respectively. The standard deviations 
provide an indication of the thickness of the shell and are likely dominated by the features on the 
western side of the shell with a wide range of distances from the center. The mean angular distances 
we derive are consistent with least-squares circle fits to each distribution independently using the 
previously mentioned common center. The least-squares fits to each transition have the added 
assumption that the distribution is circular, thus we have chosen to report only the mean distance 
from center. The computed angular distances translate to linear shell distances from center of 
4.0 ± 0.6 AU and 3.7 ± 0.6 AU at the assumed distance of 422 ± 37 pc. The standard deviations 
are likewise 0.6 ± 0.1 AU and 0.7 ± 0.1 AU for the v = 1 and v = 2 masers respectively. 

In principal, a comparison of the relative spatial locations of the v = 1 and v = 2 masers should 
allow us to comment on possible pumping mechanisms (i.e. collisional, radiative, or combination) 
as in Desmurs et al. (2000). However, because the v = 2, 42.8-GHz masers are mostly confined to 
a small region on the NW side of the shell, it is difficult to unambiguously determine the relative 
shift between the two maser shells. Although the mean distance from center we determine for the 
v = 2 masers is slightly smaller that the v = 1 mean distance, the two values are consistent to 
within the errors. 



4.3. Comparison of SiO maser distribution and photospheric stellar diameter 

Figures 4 and 5 also provide a comparison of the distribution of the SiO maser spots with our 
obtained angular size of the photospheric disk at the same epoch. The average distance of the maser 
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spots from the center of their distribution for the 43.1 GHz and 42.8 GHz transitions (rsio = 9.4 mas 
and rgio = 8.8 mas) appears at 1.7 R* and 1.6 R* respectively when compared to our estimate of the 
-fT-band UD diameter. When compared to our estimate of the continuum diameter these values are 
slightly larger at 2.0 R* and 1.9 R* respectively. The widths of the distributions are approximately 
0.3 R* for the masers in both transitions. As a reference the angular photospheric radius is listed 
in the lower panel of Figures 4 and 5 and is indicated by a colored circle at the center of the maser 
distribution. The position of the near-infrared stellar disk as indicated on the plots is assumed to 
coincide with the center of the SiO maser distribution as discussed in Sect. 4.2. The true location 
of the star relative to the masers is still unknown. We note that the above estimates are based 
on angular sizes derived for the SiO maser emission and the near-infrared photosphere, and hence 
are independent of the distance to S Ori. This result is consistent with theoretical estimates by 
Humphreys et al. (2002) based on a stellar hydrodynamic pulsation model combined with an SiO 
maser model. They obtain values between 1.7 R* and 2.1 R*, depending on the variability phase. 

Our result is consistent with a SiO maser ring radius of approximately 2 R* as determined by 
Cotton et al. (2004) for several Mira variables. Recently, however, Monnier et al. (2004) updated 
the stellar diameter estimate of Greenhill et al. (1995) for the supergiant VX Sgr and found that the 
SiO masers lie at a greater distance from the photosphere, 3.9 R*, rather than the 1.3 R* determined 
by Greenhill et al. (1995). It is striking that the stellar diameter for the other supergiant observed 
by Monnier et al. (2004), NML Cyg, also indicates a SiO maser ring distance of ~4R* when 
compared to the SiO maser shell diameter obtained by Boboltz & Marvel (2000). It is unclear 
whether these larger distances are an artifact caused by the non-contemporaneous measurements 
or, whether they indicate an inherent difference between Miras and supergiants. In earlier studies, 
measured sizes for SiO maser shells have been compared to stellar sizes based on diameters found in 
the literature (e.g. Boboltz, Diamond & Kemball 1997; Boboltz & Marvel 2000; Hollis et al. 2001). 
These diameters are often widely spaced in time and stellar phase from the measurements of the SiO 
maser shell. Sometimes the stellar diameters have not been measured at all and comparisons are 
made with even less precise diameter estimates based on the luminosity of the star (e.g. Diamond 
et al. 1994; Colomer et al. 1996; Sanchez Contreras et al. 2002). These examples demonstrate that 
it is desirable to compare SiO maser ring diameters with stellar diameter estimates, and highlight 
the need for contemporaneous multi-wavelength observations. 



5. SUMMARY 

Using the VLB A and VLTI/ VINCI, we have undertaken a coordinated multi- wavelength study 
of the Mira variable S Ori. The VLBA observations resulted in the first interferometric images of the 
v = 1, J = 1 - 0, 43.1-GHz and the v = 2, J = 1 — 0, 42.8-GHz SiO transitions. These images show 
that the masers lie in a clumpy ring-like distribution. Analysis of the VLBA images provided mean 
distances from center of 9.4 mas (4.0 AU) and 8.8 mas (3.7 AU) and standard deviations of 1.4 mas 
(0.6 AU) and 1.7 mas (0.7 AU) for the v = 1 and v = 2 SiO masers respectively at a stellar phase 
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of 0.73. From near-infrared K-band visibility measurements made with VLTI/VINCI, we derived 
UD diameters that decreased over time from ~ 10.5 mas at stellar phase 0.80 to ~10.2 mas at phase 
0.95. After extrapolating the UD diameter to the phase of our VLBA observations and considering a 
correction from UD diameter to continuum photospheric size, we determine a photospheric diameter 
of 0phot ~ 9.2 mas, corresponding to a linear radius of -Rphot ~ 1-9 AU. 

Because these observations of S Ori were closely spaced in time in a coordinated effort between 
the two instruments, we are able to relate the stellar diameter with the size of the SiO maser shell 
without the uncertainty caused by the inherent variability of the star. Thus we can conclusively 
say that the SiO masers lie relatively close to the stellar photosphere at a distance of ~2 R* at the 
time of our observations. 

The observations presented here represent the first in a series of coordinated VLBA/VLTI 
experiments to study long-period variable stars. The present VLTI observations were obtained 
during the commissioning period of the VLTI with the commissioning instrument VINCI and the 
small (40 cm) test siderostats. Future observations will utilize the scientific VLTI instruments in 
the mid-infrared (MIDI) and near-infrared (AMBER) with the 8 m Unit Telescopes and 1.8 m 
Auxiliary Telescopes. The upcoming near-infrared VLTI instrument AMBER will allow interfero- 
metric measurements in the near-infrared J, H, and K bands with a spectral resolution of up to 
10 000, and will provide closure phases. This will enable us to directly measure the the continuum 
stellar diameter and to study the conditions on the stellar surface, including possible asymmetries 
and surface inhomogeneities. 

The near-infrared results are based on public data collected at the ESO VLTI, Paranal, Chile, 
in the framework of our P70 shared risk program. We acknowledge support by the ESO DGDF. We 
are grateful to M. Scholz for valuable comments with respect to definitions of Mira star radii. We 
thank as well T. Driebe and K. Ohnaka for helpful discussions. MW is grateful for hospitality at the 
U.S. Naval Observatory. DAB is grateful for the hospitality at the European Southern Observatory. 
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Table 1. Properties of the observed VLTI/VINCI calibration stars (Borde et al. 2002). 
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Table 2. Calibrated VLTI/VINCI visibility values. 
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Table 3. VLTI/VINCI observations and derived K-band uniform disk (UD) diameters with 

associated formal errors. 



Date 




Mean 
JD 


#of 
series 


Mean 
Phase 


U UD 

(mas) 


(mas) 




Jan. 25 - 


30 


2452666.8 


17 


0.80 


10.52 


0.03 


0.50 


Feb. 04 - 


09 


2452678.8 


10 


0.83 


10.33 


0.07 


0.55 


Mar. 16 - 
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2452718.0 


12 
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9.98 


0.06 


3.91 


Mar. 27 - 


31 


2452728.1 


12 


0.95 


10.16 


0.03 


0.13 



Note. — The observation dates were grouped into 4 bins of 5 days 
each. 
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Figure Captions 

Fig. 1. — Near- infrared if-band squared visibility amplitudes of S Ori and best fitting uniform disk 
models for our four VLTI/ VINCI epochs. 

Fig. 2. — Total intensity image of the v = 1, J = 1 — (red) and v = 2, J = 1 — (blue) SiO maser 
emission toward S Ori. The images represent maximum pixel values over the LSR velocity range 
from +10.1 km s -1 to +25.4 km s -1 plotted as contours. Contour levels are 1, 2, 4, 8, 16, 32, 64 
times the 3a off source noise of 0.1 mJy beam -1 . Peak flux densitiy is 4.36 Jy beam -1 for the 
v = 1, J = 1 — line and 5.35 Jy beam -1 for the v = 2, J = 1 — transition. Synthesized beam 
sizes are 0.58 x 0.17 mas at a position angle of -15.4° for v = 1, J = 1 — and 0.52 x 0.17 mas at a 
position angle of -15.7° for v = 2, J = 1 — 0. 

Fig. 3. — Obtained S Ori i^-band uniform disk diameter as a function of date and variability phase. 
The dashed line denotes the best fitting linear functions to all of our values and to the two closest 
points to our VLBA epoch. Our VLBA and VLTI variability phases are indicated by arrows. 

Fig. 4. — LOS velocity structure of the v = 1, J = 1 — SiO maser emission toward S Ori. The 
top panel shows the spectrum formed by plotting maser intensity versus velocity, color coded in 
2 km s -1 velocity increments from redward (left) to blueward (right). The solid line in the top 
panel represents the scalar-averaged cross-power spectrum averaged over all of the VLBA antennas. 
The bottom panel plots the spatial and velocity distribution of the masers. The color of each point 
represents the corresponding velocity bin in the spectrum and the size of each point is proportional 
to the logarithm of the flux density. Errors in the positions of the features are smaller than the data 
points. The dashed circle is based on the mean angular distance of the SiO masers from the center 
of the distribution. The colored circle in the center shows the angular size of the photosphere as 
determined from our VLTI JT-band measurements. 

Fig. 5. — LOS velocity structure of the v = 2, J = 1 — SiO maser emission toward S Ori. The 
top panel shows the spectrum formed by plotting maser intensity versus velocity, color coded in 
2 km s -1 velocity increments from redward (left) to blueward (right). The solid line in the top 
panel represents the scalar-averaged cross-power spectrum averaged over all of the VLBA antennas. 
The bottom panel plots the spatial and velocity distribution of the masers. The color of each point 
represents the corresponding velocity bin in the spectrum and the size of each point is proportional 
to the logarithm of the flux density. Errors in the positions of the features are smaller than the data 
points. The dashed circle is based on the mean angular distance of the SiO masers from the center 
of the distribution. The colored circle in the center shows the angular size of the photosphere as 
determined from our VLTI K-band measurements. 
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